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Abstract 

Analyzing two simple experimental situations we show that from Newton's law 
of gravitation and Special Relativity it follows that the motion of particle in an 
external gravitational field can be described in terms of effective spatial fields 
which satisfy Maxwell-like system of equations and propagate with the speed 
of light. The description is adequate in a linear approximation in gravitational 
field and in a first order in j(? . 
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I. INTRODUCTION 



In this paper we discuss how the gravitational interaction for an objects with non-zero 
velocities can be described in terms of effective spatial fields. Namely, we will show that the 
force acting on a particle which moves in an external gravitational field is given by the expression 

F = mg + mv x B, (1) 

where g is the gravitational field accounting for Newton's gravitational law for a particle at rest 
(see below (|^)) and the effective field B, appearing due to the Special Relativity effects like the 
magnetic field in electrodynamics, satisfies the relation (below we assume that g -field is time 
independent) 

carlB^-nlt.e. fM^-rif}dA (2) 

In d^) (j) is the flow of the unit of mass per unit of time and dl and dA stand for the line and 
area elements. 

We will show that 

V = 3 

where Gn ~ 6,66 ■ 10~^^ Nim?kg~'^ is a Newton's constant and c is the speed of light. The 
small value of 47r G^v/c^ explains why, in contrast to their electromagnetic counterparts, "grav- 
imagnetic" effects caused by B are weak for moderate values of masses and velocities. 

We consider the case of a weak gravitational field and neglect higher powers of jc?. In 
terms of General Relativity we would say that the curvature of space time nearly vanishes, so 
that Special Relativity can be apphed with accuracy 1 — (700 — 1, (700 being the 00 component 
of the metric tensor g^^,- When (700 ^ 1 the curvature is almost zero, but nearly vanishing 
deviation from the flatness of space-time still leads to noticeable acceleration, described with 
accuracy jc? by Newton's gravitational law (Dirac 1976, Landau 1962). 

Below in section 11 we consider two gedanken experiments with the point particle and mass 
flow. We will evaluate t] comparing the result of the flrst experiment with the expression (|^) 
and then show that the results for the experiment 2 are described by ([l|) with the value of r] 
obtained from the analysis of experiment 1. We derive equations for the flelds g and B which 
are similar to Maxwell's equations for the electromagnetic fleld. 

In section 111 we summarize our results and discuss the hmitations of the suggested ap- 
proach. 

The concept of spatial gravitational forces modelled after the electromagnetic Lorentz force 
has a long history and many names associated with it (Holzmuller 1870, Tisserand 1872, 
Zel'manov 1956, Cataneo 1958, Bel 1959, Jantzen 1990, Damour 1991, Bonnor 1995, Bra- 
ginsky 1997, Mashoon 1997, Dunsby 1997, Maartens 1997) . In this paper we consider spatial 
gravitational flelds in the most elementary way and show that even in such a simplifled scheme 
gravitational phenomena can be analyzed with the accuracy oiv^ without invoking equa- 
tions of General Relativity. 
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II. EXPERIMENTS WITH POINT PARTICLE AND MASS FLOW AND THE 

FIELD EQUATIONS 



A. Description of Experiments 

In this section we consider two gedanken experiments: 1) point particle moving between 
two infinite pipes which carry a mass flow, and, 2) one pipe, moving towards the particle. We 
analyze these two experiments using only Special Relativity and Newton's gravitational law 

Fi2 = Gn — 3 — ri2 (4) 

The set up for experiment 1 is two infinite, straight, massless pipes on a plane and a point 
particle between them. Each of the pipes is parallel to the y axis, crossing the x axis at x = ±6, 
and each carries a fluid which flows with velocity Vf relative to the pipe. The fluid in the 
left pipe moves in the positive y direction, the fluid in the right pipe moves in the negative 
y direction. The linear density (we neglect the pipe cross- sect ion) of fluid at rest is a. The 
particle with mass m moves along y axis with the velocity Vp = Vf. 

Let us calculate the net force F acting on the particle in the frame of reference where particle 
is (momentarily) at rest - reference frame comoving with the particle. In this frame of reference 
the particle lies between two continuous mass flows with the linear densities cr and a'-fy where 
'jv is a Lorentz factor accounting the relativistic length contraction 



V 
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7y = (1 - ^)-'/' (5) 



In d^) = —j^^^Tj^ is the velocity of a fluid from a right pipe in this frame of reference. The 

forces exerted from the pipes have only an x components and the straightforward calculation 
leads to the following expression for the force acting on a particle in the frame where particle 
is at rest: 

2GNma AGNmavj 
F=^^(7.-l)-^^^ (6) 

To obtain (^ we use expression (1), e.g. the force exerted from the left flow (second term of 
(I)) is 

/■+°° GNmab ma 

The magnitude of the force exerted from the right flow is given by (|^) using substitution 
a — *• ajv (since we consider inflnite pipes there are no any boundary effects caused by the 
endpoints of fluid) and the flnal result is the expression (|^). 

When the particle is at rest relative to the flow the net force is zero, as it follows from 
(^). Below we will show that B ^ Vf (see (PUf)) and that expression (^ is reproduced by the 
w-dependent term of the equation (|l|). 

In the second experiment, the particle with coordinates {x,y,z) is at rest and the pipe 
oriented along y axis moves towards the particle with the velocity v = (— fp,0,0). The fluid 
with density p moves with the velocity Vf = {0, —vj, 0) relative to the pipe. 
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We omit lengthy but straightforward calculations and report only the results for a - the 
acceleration of the particle: 



_ dvy 2GNO'VpVfX 

where x is times the distance from the pipe to the particle in the reference frame where 
particle is at rest. Expressions are obtained by integration similar to that ([^ which is based 
on (^), taking into account the length contraction of a small element of fluid with mass pdy 
and neglecting higher orders of v'^/c?. 



B. Evaluation of r] from the results of Experiment 1 



To obtain the value of r/, appearing in the relation (^) it is necessary to consider the problem 
in a reference frame where particle has a non-zero velocity. The simplest solution is provided 
by the original frame of reference described previously when particle and the fluid from a left 
pipe move along y axis with velocities Vp and Vf correspondingly, Vp = Vf, and the fluid from 
the right pipe moves with the velocity —vj. 

From the symmetry arguments it follows that the net field g is zero (particle is between 
two sources with the same linear densities (J7), so only the second term of (|I|) contributes. 
According to (|l|) the force acting on the particle is directed along x axis and its magnitude is 
mVpB. This expression already assumes that Bvf = Bvp = 0, i.e. B is perpendicular to the 
pipes just as the magnetic field given by Ampere's circuital law is perpendicular to the current 
(Landau 1962). In experiment 1 we have dg/dt = (pipes are at rest), so we can use (0) to 
calculate the value of B. Particle is at rest relative to the fluid in the left pipe, so it "feels" 
field B generated only by the pipe from the right. Using the integral relation (Q) and for the 
relative velocity V = —j^^j^ we obtain 

2vf 

2TxhB = —r]a'jV = -qa — ^ 2r]avf, (9) 

1 + ^ 

i.e. 



B ^ riavf/irb (10) 

and the force acting on a moving particle in the original reference frame will be (taking into 
account Vp = Vj): 

rjmavj 

= mVpB = (11) 

no 
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In order to establish the value of rj we have to compare F and JF - force calculated in a two 
reference frames. Reference frames move relative each other along y axis with velocity V and 
the force is directed along x axis. Thus we need Lorentz transformations in its vectorial form: 

, = .(t - (^); ? = 4r'f- VtUl- (12) 



From (O) and the expression for the force 



it follows that up to order v'^ jc^ (let us remind that F itself is of order f^/c^) we have F = T . 
Therefore, equating expressions for F and T we obtain 

'^''^'l^o(')\ (14) 



C. description of Experiment 2 in terns of g and B 

After the value of ri is established we demonstrate that the results for experiment 2 are 
described by (|1]). To do so, instead of straightforward but lengthy arguments we assume that 
when dg/dt ^ (in the reference frame used to describe experiment 2, the pipe has non-zero 
velocity, so g is now time dependent) expression (^ is modified as 

-* 1 dg -> f ^ ^ f 1 (^9 ^ ^ 
curlB = — — - T](j), i.e. j Bdl = J (^"^ - V(p)d^ (15) 

We will verify that (HD and ( |l5l) lead to the same expressions for the force as did calculation 
using only Newton's law and Special Relativity. 

The components of B can be calculated from (plSf). For g we use g = a with a given by (|^) 
where dg/dt = (dg/dx) ■ (dx/dt), and assume that fields vanish at infinity. The choice g = a 
is justified since in the original frame of reference particle is at rest and its (instantaneous) 
acceleration is defined by a term "of (|I]). Using rj = AttGn/c^, in a leading approximation 
in v'^/c^ we obtain: 



B 



2GNcrVfZ 

C^{X'^ _|_ ~;2^ 

2G]\faVpZ 



- -4^f^ (16) 



Now we are in a position to check that experiment 2 can be described by F = mg + mv x B 
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When the particle is at rest in the reference frame used to describe experiment 2, F = mg, 
which is trivially consistent with the equations since we have defined g = a. 

To account the effect caused by the "B - term" of (0), we consider the case when particle in 
the original frame of reference moves with the non-zero velocity u = {u, 0, 0) . The y component 
of (|l]) is: Fy = mgy — mBzU. Direct substitution for Bz (and gy = ay) results in 



To compare (0) with the expression calculated in the framework of Newtonian approximation, 
Fy = may, we need the value of ay. From (P), acceleration in case of = 0, it follows that 
when a particle has non-zero u = {u, 0, 0), the value of ay can be obtained by substitution of 
Vp + u for Vp in (H) : Vp ^ Vp + u which leads to 

2GnO'(vp + u)vfX 

Fy = may{vp ^vp + u)= m ^^^^^ ^ (18) 

As it is clear, (0) and (p!8D agree. 

Next we consider the z component. Using (|16D for By we obtain: 



mgz + muBy = -m——-{\ + + - m J u (19) 



^2Gjv^ 2Gj^aVpZ 
'x^ + z^^ 2c2 2c2^ ^c^{x^ + z^ 

Since the particle moves with the velocity u = {u, 0, 0) we have a^ = Fz/'jufn (see (|13[)): 

2GNcrz vl vj 2G^aVpZ 



a;2 + ^2v 2c2 2&' cHx^ + z^] 



where 7"^ = y 1 — u'^/c^ 1 — u'^/2c^. 

Now we have to compare this expression with the one for from (^, calculated from 
Newton's law and Special Relativity - acceleration of a particle in the reference frame comoving 
with the particle . We replace Vp Vp + u in the expression (^) for a^ to obtain: 

2GNaz. vl vl v„u -u^ 

Expression ( pi]) gives the acceleration of the particle in a reference frame moving along the x 
axes with velocity Vp + u relative to the pipe. On the other hand, expression (|20| ) corresponds 
to the acceleration in the reference frame moving along the x axes with velocity Vp relative to 
the pipe. To compare (^) and (^), we use (|T^ to transform acceleration from the reference 
frame used in ( pTf ) to that used in (pOl): — > a'^'j~'^. This transformation introduces the term 
~u^/c^ so that the two expressions for acceleration now agree. 

A similar analysis for the case when the velocity of the particle is along the y axis, u = 
(0, u, 0) again confirms that the expression F = mg + mv x B can be used to describe the 
motion of a particle in a gravitational field. 
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D. Equations for g and B 

Besides equation ([T5|) which was postulated (and subsequently verified to describe experi- 
ment 2 in a self-consistent way) it is possible work out two more relations for the fields B and 
9- 

First of all, from the expressions (|TB|) it follows that dBx/dx + dBy/dy + dBz/dz = 0, i.e. 



divB = (22) 

— * 

Next we compare curl g and dB/dt. For the y component we obtain: 

dBy d f 2GNpvpZ \ 2GNpVp2zx 



dt dt \ c^{x^ + z^) ) c^{x'^ + z'^Y 



(23) 



Straightforward calculation of the y component of curl g (as before, we take g = a, for a see 
(S)) results in 



dgx dg, _ d / 2GnPx . ^} ^\ d ( 2GnPz \ _ 2Gnp2zx 

dz dx dz\ {x^ + z^y 2c^') dx\ {x^ + z'^))'^ {x^ + z^f d^' ^ ' 

i.e. {curl g)y = —dBy/dt. 

Consideration of the x and z components show that the relation 

dB 

curl q = ——— (25) 
dt ^ ^ 

is valid. Also, from the definition of g we have div g = AttGnp where p is a regular three 
dimensional density. 

Summarizing, the equations for g and B are as follows: 

div g = AuGnp, curl g = —dB/dt 

div B = 0, curl B=\^4- (26) 

ot c^ 

where p is the mass density and j is the mass density fiow. In case of experiments 1 and 2 

j=(P = pVf. 

It is now straightforward to obtain the wave equations for the case p = 0, j = 0: 

1 d^g 2^ 1 d^B 2 

i.e. free waves propagate with speed of light. 
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III. DISCUSSION 



We have demonstrated that the gravitational force acting on a point particle with mass m 
and velocity v is given by the expression 

F = mg + mv x B (28) 

— * 

with g and B satisfying the system of equations similar to the Maxwell equations. 

The approximation we used is that gravitational field is weak enough so that space-time 
is approximately euclidean and the velocities are small enough so that higher powers of v"^ /(? 
are negligible. In the framework of this approximation the force obtained from Newton's law 
(I) and the Special relativity is described by (^81) , i.e. motion of particle is given by an 



expression similar to the Lorentz force for a charged particle in an external electromagnetic 
field. The similarity is caused by neglecting the effects of self-interaction for gravitational field, 
corresponding to a non linearity of Einstein's equations. In case of classical electromagnetism 
the linear approximation to field equations is well justified in a sense that phenomena with 
characteristic action substantially exceeding h being the Planck's constant, are described by 
Maxwell's and Lorentz's equations (Landau 1962) and in electromagnetic phenomena quantum 
effects manifest themselves earlier than effects caused by a non linear corrections to Maxwell's 
equations. Intuitively it becomes clear when one compares electron's Compton wave length 
Vq = h/mc and its classical electromagnetic radius r^. = e^/mc^: from the value of the ratio 
fq/Tf, = hc/e^ ~ 137 it follows that the quantum effects, namely the pair production occurs at 
a distance which is 137 times more than the distance at which the classical field singularities 
become relevant and when it becomes necessary to modify classical theory, e.g. to introduce 
non linear terms in field equations. 

Theory of gravity provides us with an opposite feature - "classical radius" = IGj^mjc^ 
appearing in the Schwarzschild's metric (Dirac 1976, Landau 1962) greatly exceeds Compton 
wavelength - r^/r^ = 2m? /Mpi, where the Planck mass Mpi ^ IQ^^gr. Therefore in describing 
the motion of bodies with m ^ Mpi it is vital to consider the exact classical equations of 
motion (Einstein's non linear equations) - quantum effects are negligible at this scale. The self- 
interaction plays a decisive role in describing basic phenomena of light deflection or precession 
of perihelion of planetary motion (Dirac 1976, Landau 1962). This features of a motion in a 
static gravitational field can not be described by (^) - the "i3 term" is absent for a static source. 
Therefore, fields g and B, describing gravitational interaction of a moving objects in a linear 
approximation, can be treated only as an effective fields and the limitation of our approach 
manifests itself in a degrees of freedom: 6 components of g and B of course are not enough to 
describe the degrees of freedom of a gravitational field. 

Despite that the relevancy of the linear approximation is questionable, approximation (^) 
can be still useful for describing interaction of moving bodies : from the General Relativity it 
follows that the exact expression for the force exerted on point particle moving in an external 
stationary field is given by expression similar to (|28|) (Landau 1962): 



F = —mc \/ln \^—goo + "mc^Z—gooV x curl Q, (29) 

where Qa = —Qai)/ Qmi ot stands for a spatial part of metric and g^^ is a metric tensor. When 
^fQo = ~1 ~ 2$/c^, where $ is a scalar potential, $/c^ -C 1, the first term of the r.h.s. of (p9|) 
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is the same as the first term of the r.h.s. of (pSl) . To reproduce the second term of (pH]) which 



includes the field B, that is to express B in terms of it would be necessary to solve Einstein's 
equations. At the moment, we know of no solutions for the Einstein's equations for experiments 
1 and 2, but based on our phenomenological consideration we believe that the equation similar 
to ( p8D can be derived from the equations of General Relativity. 

Let us note that though fields g and B satisfy wave equations (|27|) , they do not transform 
as an antisymmetric tensor of rank 2, i.e. they do not transform as the electromagnetic field 
Ffj^i, ~ {E, H). If one attempts to postulate that the exact expression for the force acting on a 
test particle is given by (^) or ( P^D then it turns out that in order to maintain expression (^) 
fields g and B transform like non-tensor quantities (Dirac 1976). This is the price one has to 
pay when attempting to describe gravitational interaction in terms of 6 degrees of freedom. The 



non-tensor feature of transformation is most transparent from (|29|): identifying g with the first 
term of r.h.s. of ( PU] ) it follows that at ^ + ^^{x) in the expression for the transformed 
g there arises extra term 

In ^/-goo, (30) 




which can not be compensated by the transformation of a "i3-term" of (PP|). Therefore eq. (^) 
can not hold in any reference frame, for any velocities. In Einstein's equations extra terms 
similar to ( PD| ) are compensated by coordinate transformations of General Relativity and as 
a result, equations of gravitational field and the requirement of general covariance form a self 
consistent mathematical scheme (Dirac 1976, Landau 1962). 

In approximation used in this paper (linearized equations and lowest order in v'^/c^) fields 
g and B transform as E and H. This statement is true only in lowest order in v'^/c^. Straight- 
forward calculation shows that (^ and (|T6|) transform exactly as E and H transform in the 
lowest order in v'^/c^. Therefore, in the framework of approximation used the description based 
on ( p8D and (|26| ) is self consistent. 

As we already have mentioned, g and B are effective fields, even from the point of view of 
classical theory. Nevertheless, equation ( pH]) can be applied to a rather wide class of phenomena 
in problem of describing the motion in external gravitational field after the fields g and B are 
known. The advantage of using (^) and (|26| ) is in their simplicity in comparison with the 
problem of solving equations of General Relativity. 
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